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Abstract 

How nanoscale uniaxial magnets respond to an alternating field is studied by 

direct numerical calculation. A nontrivial oscillation of the magnetization is 

found, which is analyzed in terms of the non-adiabatic transition due to the 

time dependent field. A new method to estimate the tunneling gap of the 

magnet is proposed. 
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The quantum dynamical behavior of nanoscale magnets have attracted interests both the- 
oretically and experimentally. [|9|-|T3H In particular, the step-like magnetization process 
of Mni2-Ac is a peculiar phenomena, due to quantum dynamics, for which various expla- 
nations have been proposed. We have also proposed a mechanism in terms of successive 
Landau-Zener-Stuckelberg(LZS) transition. |T3|-|T5|] We pointed out that the non-adiabatic 
transition between nearly degenerate levels gives essential the mechanism for the steps and 
that the velocity of the changing field, c, is an important parameter because the transition 
probability is very sensitive to c. 

For the LZS transition, so far mainly a linear time dependence of the field, H(t) = 
ct, has been considered, although an effect of oscillating field has been studied in a very 



different context. JT7J In the present letter, we study the quantum mechanical response to 
an alternating field for a simple model of a uniaxial magnet. Here we study the transverse- 
field Ising model: 

n = -J £ <^ -rX>f- H(t) J2 (i) 

<i,j> i=l i=l 

where 

H(t) = H cos cut. (2) 

Throughout this letter we take J as a unit of energy and put it unity. Here we show only 
results for a system of four spins (L = 4) subject to periodic boundary conditions and 
T = 0.5. For other choices of model parameters, we find qualitatively similar behavior. In 
Fig. 1 the energy levels of the model are shown as a function of the field H = H(t). Only 
the 8 lowest states, some of which are degenerate, are shown. When the energy gap at 
H — is small, the lowest two levels are located far below the other levels. When we take 
the initial state to be the ground state, the system can be regarded as a two level system as 
far as Hq is small and does not cause the second scattering to the higher levels. Successive 
non-adiabatic transitions to higher levels have been studied in Ref. 0. 

If H Q is very small, we may use the Kubo formula to study the linear response, 



(M z {t)) = Re( X (u)H e 



iujt\ 



(3) 



where xi 00 ) is the dynamical susceptibility 

X(u) = - lim - / Tr[M z ,M z (t)]e- m °- lujt - et dt/Z, 

e^+o h Jo 

which is given at zero temperature, namely (3 — > oo, as 

2hu 



(4) 



xM = £|<g|m, 



v- 



— iir5(Ee — Eq — hu) 



(5) 



(E e -E G )*-(hw) 2 

where 7i is H. with H = 0, M z = Y^i a t an d ^ = Tre^ 137 " 10 . E G and i?^ are the energy 
of the ground state and the £-th excited state, respectively. In this formula the Zeeman 
term H(t)M z is treated as a perturbation and relevant frequencies are only those due to 
the energy gaps at H = 0. On the other hand, in the present paper we are interested in 
the phenomena due to the non-adiabatic transition where H(t)M z can not be treated as a 
perturbation. Thus even when we say H is small, H Q is still 0(1) and not small enough to 
be treated clS cL perturbation. 

The probability for staying in the ground state when the field changes the sign is given 
as 

/ n(AE) 2 ^ 



p = 1 — exp 



where c is the velocity which is given by 



4c|M 



(6) 



HqU) 



(7) 



H(t)=0 

and M is the ground state magnetization near H = 0, 



M = lim M(H). 



(8) 



For T = 0.5, the energy gap at if = between the ground state and the first excited state, 
AE, is 0.03549 and |M | ~ L = 4. 

The time evolution of system is given by 



\t) = e - l Io n(s)ds/h \0), 



(9) 



where |0) is an initial state which is chosen to be the ground state of the model for H = H(0) 
and the exponential denotes the time-ordered exponential. We solve Eq.(|9|) making use of the 
4-th order decomposition proposed by Suzuki. ||18|,|19j Hereafter we put h = 1 for simplicity. 

As has been shown in the previous studies, Eq.(|6|) is confirmed by the simulation results. 
From Eq.([jp, p = 0.0062 for uo = 0.2 and H = 0.2. In the simulation we calculate the 
overlap between the ground state and the first excited state 

x(t) = \(G(t)\t)\\ (10) 

where \G(t)) is the ground state for H = H{t). After a half period, t = ir/ui, we find 
x(t = 7r/ou) = \(G(t)\t)\ 2 = 0.0063, which confirms the LZS prediction. 
In Fig. 2 we show the time dependence of the magnetization, 

"»(*) = ME °?l*>, (ii) 

i 

and observe a gradual relaxation due to the successive non-adiabatic transitions. When we 
continue the simulation, a sinusoidal motion is found as shown in Fig. 3, 

m(t) ~ cos(nt), (12) 

where x(t) and H(t) are also shown. The period of this sinusoidal motion does not correspond 
to an eigenfrequency of the system nor to the period of the external field. Actually when 
we change the amplitude of the field H Q the period of the magnetization changes as shown 
in Fig. 4(a). The dependence of the period on u is also shown in Fig. 4(b). 

Although the dependence of Q on H in Fig. 4 seems irregular, we find a rather regular 
dependence when we plot the frequency Q as a function of H , shown in Fig. 5. 

Let us study the dependence of Q on Hq and uj. The time-evolution during the field 
changes from H to — H is given by 




which is 2x2 unitary matrix as far as H is small and only the lowest two states takes 
dominant role. Here we take the ground state \G) and the first excited state |1) for the 
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initial state as the basis. Because the probability p is known as Eq.(^|), we take the following 
form of X 

X=[ t , (14) 

where 9 and are unknown phases which depend on H Q and u. After t = ir/u, an inverse 
process is taken. If we change the sign of ^-component of the spinsfj the time-evolution 
during the field changes from — Hq to Hq is identical to X, because the motion of the 
Hamiltonian is identical. Thus we only have to change the basis. The ground state \G') and 
the first excited state for t — ir/u is generally expressed as a linear combination of \G) 
and |1): \G') = a\G) + b\l) and = c\G) + d\l). Let the transformation matrix be Q. 
Thus the second half time-evolution X' is expressed as 

X' = exp ( -z H(s)ds) , = Q~ X XQ. (15) 

When the scattering region of H(t) is very narrow, which is the present case as shown in 
Fig. 1, we may take 

«=(;;). d«) 

Combining Eqs. ([l4|) ,([T5D and (0), the time-evolution operator for one period is given as 



C = Q- X XQX 
e 2ie p + (i _ p ) e # ( e ^+i0 _ e ~ie+2i^ Jp{i- p ) 



(e ie - e- ie+l<t, )^jp(l - p), e~ 2ie p + (1 - p)e** 

The eigenvalues, A±, of C are given by 



(17) 



\± = (q±ijl-<?)e i + (18) 

where 



For this change we can use the unitary transformation (a x , a y ,a z ) — ► (a x , —cr y ,—a z ) for all sites. 
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q = 1 — p + p cos(a) , a = 29 



(19) 



Now we put 



(20) 



tan 



where 2tt/uj is the period of the field and we take Q to give the frequency in Eq. (|T2|) . Here 

(21) 

V u; / q 1 — p(l — COS a) 

When p 1 , 

ttO / 

(22) 



7iQ\ _ y/1 - q 2 _ V 2 ^ 1 ~~ c °s«) -P 2 (l - cosa) 
lo J q 



2p(l — cosa). 



The evolution of the wave function is generally given by 



i 



c+e 



y+ 



X _ 



+ c_ 



+ c„e 



( 










) 


y- 





(23) 



where t (x±,y±) are the eigenvectors and c± are coefficients. Thus the probability being in 
the ground state after n periods, x(2nn/u), is given by 



x(27in/uj) = \c + x + e l * n + c_x_e" 
= a + ocos( n + 7). 



(24) 



When p <C 1, x + ~ a;_ ~ y + ~ — y_ ~ l/\/2 and c + ~ c_, a ~ 6 ~ 1/2 and 7 1, which 
explains the time-evolutions in Fig. 4. 

Here let us estimate Q from x(^) and x(— ). From Eq . (|T4"1) and Eq . (|P7D , 



and 



x(-) = p, 
a; 



2ir 

x(-) = 1 - 2p(l - cosa) + OQ» 2 ). 

UJ 



(25) 



(26) 



Thus 



R 




2(1 — cos a). 



(27) 



In Fig. 6 we show R estimated from numerical calculations of x(ir/u) and x(2h/oj) for various 
values of Hq, which agrees with the Hq dependence of Q shown in Fig. 5 (Q = \fpR x uj/2tt). 
In Fig. 5, fi max = yjp x uj/tt is shown by a dotted line. We find that y/p x uj/tt gives the 
envelope of tt(Ho) and confirm that the relation Eq. (|22"D holds. 

Generally we know neither a nor p. Even in such situation, we can estimate fi max by 
observing Q for various Hq and u. Alternately, from f2 max , we can estimate p by 



and therefore also DeltaE by making use of the relation Eq.^). The present analysis is 
good for any periodic function of H(t), not necessarily cos(ut). We have confirmed the same 
behavior for a piecewise linear function (i.e., zigzag function) of H(t). 

The present oscillation of M(t) is due to the non-adiabatic transition and is a peculiar 
property of quantum dynamics with time dependent field. The present mechanism is so 
simple that it would be applicable for many cases where the non-adiabatic transition takes 
place and that we hope such nontrivial oscillation would be observed in an experiment of 
nanoscale systems. 

The present study is partially supported by the Grant-in- Aid for Scientific Research from 
the Ministry of Education, Science and Culture. 
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FIGURES 

FIG. 1. Fig.l Energy levels for 7 = 0.5. 

FIG. 2. Fig.2 Time dependence of magnetization, M[t) is shown by a solid line. H(t) + 2.5 is 
shown by a dotted line. 

FIG. 3. Fig.3 Nontrivial oscillation of magnetization M(t). x(t) and H(t) are also shown. 

FIG. 4. (a) Time dependences of magnetization for various amplitudes, Hq. (b) Time depen- 
dences of magnetization for various frequencies u. 

FIG. 5. Hq dependence of Q. £l max is shown by a dotted line. 

FIG. 6. Hq dependence of R. 
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M(t) r=0.5, a>=0.2, H =0.2 




Fig. 2 



M(t) r=0.5, a>=0.2, H =0.2 




Fig. 3 



Resonance o>=0.2, H =0.1 8,. .0.26 
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Resonance o»=0.1 , 0.2, ..0.5, H =0.20 




Fig. 4b 
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